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Abstract 

We study bulk fermion fields in various multi-brane models with localized gravity. The 
chiral zero mode that these models support can be identified as a right-handed sterile neu- 
trino. In this case small neutrino Dirac masses can naturally appear due to a localization of 
the bulk fermion zero mode wavefunction, in an analogous way to graviton, without invok- 
ing the see-saw mechanism. The conditions and the options for localization are discussed 
in detail. It is shown that, considering a well motivated five dimensional mass term, the 
localization behaviour of this mode can resemble the graviton's at least in a region of the 

parameter space. As a result, the " H h", " + +" models can support, in addition to the 

ultralight graviton KK state, an ultralight localized and strongly coupled bulk fermion KK 

mode. We find that there are severe constraints on the parameter space of " H V" and 

" + +" models if the neutrino properties resulting from this light fermion state are to be 
reasonable. Furthermore, in the case that also the Bigravity scenario is realized the above 
special KK mode can induce too large mixing between the neutrino and the KK tower 
sterile modes restricting even more the allowed parameter space. 



1 s. mouslopoulos@physics.ox.ac.uk 



1 Introduction 



The study of bulk fermion fields, although not something new 0, turns out to be of 
particular interest in the context of brane-world scenarios both in the case of models with 
large extra dimensions |3], |], |5| (factorizable geometry) and in models of localized gravity 
H 0, |8| (non factorizable geometry) since they can provide possible new ways to explain 
the smallness of the neutrino masses, neutrino oscillations and the pattern of fermion mass 
hierarchy. 

In the context of string and M-theory, bulk fermions arise as superpartners of gravi- 
tational moduli, such as, those setting the radii of internal spaces. Given this origin, the 
existence of bulk fermions is unavoidable in any supersymmetric string compactification 
and represents a quite generic feature of string theory^. This constitutes the most likely 
origin of such particles within a fundamental theory and, at the same time, provides the 
basis to study brane-world neutrino physics. 

In the traditional approach the small neutrino masses are a result of the see-saw mech- 
anism, in which a large right-handed Majorana mass Mr suppresses the eigenvalues of the 

neutrino mass matrix leading to the light neutrino mass m v ~ ^y 7 " " ■ The neutrino mix- 
ing explanations of the atmospheric and solar neutrino anomalies require that Mr to be a 
superheavy mass scale > 10 10 GeV. 

In the case of large extra dimensions || [|, || , despite the absence of a high scale like 

Mr (since in such models the fundamental scale can be as low as 1 TeV), small neutrino 

masses |TT], [II]] (Dirac or Majorana) can arise from an intrinsically higher-dimensional 

mechanism. The idea is that any fermionic state that propagates in the bulk, being a 

Standard Model (SM) singlet can be identified with a sterile neutrino which through it's 

coupling to the SM left-handed neutrino can generate small neutrino mass. In the case 

of factorizable geometry, the smallness of the induced masses is due to the fact that the 

coupling is suppressed by the large volume of the internal bulk manifold. In other words, 

the interaction probability between the bulk fermion zero mode, the Higgs and Lepton 

doublet fields (which are confined to a brane) is small because of the large volume of bulk 

compared to the thin wall where the SM states are confined, resulting a highly suppressed 

coupling. In the context of these models one can attempt to explain the atmospheric and 

solar neutrino anomalies (see e.g. |TD|, O [13], |TJ], 

2 However, note that brane-world models with non factorizable geometry have not yet been shown to 
have string realizations. For string realizations of models with large extra dimensions see Ref.[|) 
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Figure 1: The right-handed fermion zero mode (solid line), first (dashed line) and second 

(dotted line) KK states wavefunctions in the symmetric " H V" model. The same pattern 

can occur for the corresponding wavefunctions in the " + +" model. The wavefunctions of 
the zero and the first KK mode are localized on the positive tension branes. Their absolute 
value differ only in the central region where they are both suppressed resulting to a very 
light first KK state. 



In the context of brane world models with localized gravity || |7|, || (non factorizable 
geometry) small neutrino masses can again be achieved, without invoking the see-saw mech- 
anism. In more detail, in this case the mechanism generating the small coupling between 
the Lepton doublet and the Higgs which live on the brane and the right-handed sterile 
neutrino zero mode is not due to the compactification volume (which is now small) but 
due to the fact that the sterile neutrino wavefunction can be localized [16-28] on a dis- 
tant brane. One may thus arrange that the overlap between this mode and the SM brane 
is sufficiently small. In this case the AdS§ geometry localizes the fermion zero mode on 
negative tension branes. Localization can occur on positive branes if a mass term of the 
appropriate form is added to counterbalance the effect of the AdS$ geometry, by applying 
the ideas presented in Ref. |T|, Such a mass term appears naturally in the case that, the 
branes arise as the limiting cases of domain walls that are created from a five- dimensional 



scalar field with an nontrivial ground state (kink or multi-kink) |2^, 24 1. In this case the 
scalar field naturally couples to the bulk fermion through an non trivial "mass" term and 
will naturally induce localization to the bulk modes on positive or negative branes - which 
depends on the vacuum expectation value of the scalar field. 

In this paper we analyze the localization of a bulk fermion the mass spectrum and the 
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coupling between SM neutrino and bulk states in the context of multi-brane worlds (e.g. 
see Ref. 0). We discuss in detail the conditions and the options for the localization in 
relation to the form of the bulk mass term. We also discuss the possibility of generating 

small neutrino masses in the context of " + H — " H h"|j, " + +" models. We determine 

for which regions of the parameter space lead to a solution of the hierarchy problem and 

generation of small neutrino masses. The study of " H h" and " + +" models reveals the 

possibility of an ultralight KK state of the bulk fermion analogous to the KK graviton in 



the gravitational sector |29| , |32| , |33| , p4fl . This is due to the fact that, in this region the 
wavefunction of the right-handed bulk fermion states obeys a similar equation to that of 
the graviton. The above fermion state, when exists, imposes even more severe constraints 
on the parameter space of these models. 

The organization of the paper is as follows: In the next section we review the general 
formalism, and discuss the form of the bulk neutrino mass term which turns out to be 
critical in determining the neutrino properties. In section 3 we review the bulk fermion 
properties in the Randall-Sundrum (RS) model @]. In section 4 we extend this to the case 
of three brane " + H — " model |33| . In sections 5 and 6 we study in detail the "Bigravity" 



three brane models " -\ V" [32] and " + +" ||34|| . In section 7 we discuss the implications 

of the realization of the "Bigravity" in relation to the presence of bulk spinors. The overall 
implications and conclusions are presented in section 8. 

2 General Framework 



Following the framework introduced in Ref.[T7]], we consider a spinor \1/ in a five dimensional 
AdS§ space-time, where the extra dimension is compact and has the geometry of an orbifold 
S 1 /Z 2 . The AdSs background geometry is described by 0: 

ds 2 = e- 2a{y) Vflu dx^dx u -dy 2 (1) 

where the warp factor cr(y) depends on the details of the model considered. For the moment 

we assume that we have a model with a number of positive and negative tension flat branes 

(the sum of the brane tensions should be zero if one wants flat four dimensional space 

3 We consider the " H h" configuration as a toy-model ignoring the phenomenological difficulties as- 
sociated with the presence of a negative tension brane p8[ since it's phenomenology is very similar to the 
" + +" which includes only positive tension branes. 

4 We will assume that the background metric is not modified by the presence of the bulk fermion, that 
is, we will neglect the back-reaction on the metric from the bulk fields. 
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on the branes) and that this function is known ( it can be found by looking the system 
gravitationally) . 

The action for a Dirac Spinor of a mass m in such a background is given by []: 



S = Jd*xJ dyVG{E* (dl - S) * + ^f*{l a , ° bc }V 



m<£«} (2) 



where G = det{GAB) — e 8ct ^. The four dimensional representation of the Dirac matrices 
in five-dimensional flat space is chosen to be: 7° = (7^ , 175) . The inverse vielbein is given 
by E£ = diag(e a ( y \ e a<yV \ e a ^ y \ e a<yV \ 1). Due to the fact that the vielbein is symmetric, the 
contraction of lOhcA Q with the corresponding term in the action gives vanishing contribution. 
The mass term is assumed to be generated by a Yukawa coupling with a scalar field $ which 
has a nontrivial stable vacuum $(y). 

A few comments on presence of the mass term are in order. The motivation for intro- 
ducing a mass term of this form comes from the need to localize fermion zero modes in the 
extra dimension. This is discussed in more details in the next paragraph. Note that the 
"localization" of the wavefunction of a state does not necessarily reflect the actual localiza- 
tion of the state, since in one has to take in account the nontrivial geometry of the extra 
dimension - something that is done when we calculate physical quantities. Also when we 
note that a state is localized on a brane we mean that this holds irrespectively if the space 
is compact or not (thus the state should be normalizable even in the non-compact case). 
In order to have a localized state (zero mode) on a positive brane it is necessary to have 



an appropriate bulk mass term. This is because |16|, |18j the background AdS$ geometry 
itself has the opposite effect favouring localization on negative tension branes. As we will 
discuss, this leads to a critical mass m cr below which the localization is still on the negative 
brane (if we consider a configuration of a positive and one negative brane), for m = m cr 
the is no localization and for m > m cr the zero mode is localized on the positive tension 
brane []. 

5 We do not include a Majorana mass term, <J' T C\I/, which is forbidden if the bulk fermion has a conserved 
lepton number. 

6 where w Ma6 = \{d p e hv - d v e hll )e a v - \{d p e av - d v e afl )e b v - \e a p e b a (d p e ca - d a e cp )e c fJ , 
7 One may ask why one should localize the fermion zero mode on a positive tension brane and not on 
a negative ? In the case of RS model it is obvious that if we demand to solve the hierarchy problem and 
in the same time to create small neutrino masses through this mechanism, one should localize the fermion 
zero mode on the positive brane. In the case of multi-brane models this is not a necessity and thus, in 
principle, another possibility (of course in this case the mass term - if needed - should have different form 
(e.g. for m — » —m)). 
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2.1 The mass term 



Let us return to the specific form of the mass term and it's generality since this will become 
important in later discussions. Since the mass term is a key element for the localization 
let us, for a moment assume that it has the form m&(y)tyty where $(y) is the vev of a 
scalar field which has a nontrivial stable vacuum (its vev does not depend on the remaining 
spatial dimensions). Now if we do a simple calculation, in the case of flat extra dimension 
(the general arguments will apply also in our background eq.(l), taking in account also the 
effects of the AdS$ geometry), we find that the above configuration implies that the zero 
mode satisfies a Schrondiger equation with potential of the form V(y) = & 2 (y) — $'{y)- In 
order to localize one state the profile of $>(y) should be such that it creates a potential well. 
The way to do this is to assume that the ground state of the scalar field has a kink or a 
multi-kink profile Jl], |2|, |23|. Although the details of these profiles depend on the form of 
the potential of the scalar field if we demand strong localization of the states, the kink 
profiles tend to 9- functions. This implies that the function $ in eq.(2) can be considered 
as an arbitrary combination of functions (compatible with the symmetries of the action). 
However, as shown in Ref.|23|| the same field $ can be used in order to create the branes 
themselves. This restricts the possible form of the mass term. If we assume that the same 
field $ creates the branes and localizes the fermion zero mode the mass term should have a 
(multi-)kink form, with $(?/) = up to a sign.0 The previous argument also supports the 
^-function form of the mass term (and not for example a tanh(y) profile) since we assume 
that the branes are infinitely thin. Note thought that in the " + +" model due to the AdS^ 
geometry on the branes the, cr'(y) function does not have just a ^-function form but it also 
involves kink profiles. 

If one assumes that the mass is generated by coupling to a scalar field, different from 
the one that creates the branes, it can have any form allowed by the dynamics and the 
symmetries of the action. For example it can take the form: — m (0(y) — 9(—y))). In this 
case, it will tend to induce localization on the brane siting at the origin of the orbifold. 
Nevertheless this will not be satisfactory option in multi-brane models where the desired 
M pi /Mew hierarchy is not generated between the two first branes (e.g. " H h" or even in 



8 Note thought, that if we choose the opposite sign for the mass term (i.e. m — * — m) the localization 
of the fermion zero mode will always be on the negative tension branes. Taking in account the fact that 
the latter occurs also, in a region of the parameter space, in the case that our mass term choice is the one 
that appears in eq.(2), we will not consider this possibility separately since we can easily, as we will see, 
generalize our results for m < (the presence of such a mass term sharpens the localization of the states 
on negative branes). 
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the " + +" model) since in these cases the will be not generally be possible to simultaneously 
solve the hierarchy problem and generate small neutrino masses (it would work though in 
the " + H — " model - or it could work in the cases where the background induced localization 
dominates but in that case the mass term would be unnecessary anyway). In any case since 
we are interested in the most economic, in terms of parameters and fields, models we will 
not consider these possibilities. 

The geometry we consider has a Z 2 symmetry ( y — > — y ). Under this the fermion 
parity is defined as: $(-y) = 75^(y) (i-e. = -^l(u), ^r{-u) = ^r(u) ) and 

changes the sign of a Lagrangian mass term of the form: m\I/\E' = iji^l^r + ^r^l)- The 
full mass term however is invariant under the Z 2 since the function cr'(y) is also odd under 
the reflections y — > — y . With this definition of parity one of the wavefunctions will be 
symmetric and the other antisymmetric with respect to the center of the orbiford . Note 
that this implies that the odd wavefunction will be zero at the orbifold fixed points (i.e. 
zero coupling to fields confined to that points). Since we would like in what follows to use 
the right-handed component in order to give mass to SM neutrinos, which could be confined 
on a brane at an orbifold fixed point, we choose the right-handed wavefunction to be even 
(i.e. non-vanishing coupling) and the left-handed to be odd. 

2.2 The KK decomposition 

It is convenient to write the action in terms of the fields: and where ^r,l = 
|(1 ± 75)^ and ^ = ^ R + ^ L . The action becomes: 




writing and \l/ l in the form: 



(y) 



(4) 



n 



the action can be brought in the form 



S = S / ^xi^ix^d^nix) - m n ip n (x)ip n (x)} 



(5) 



n 
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provided the wavefunctions obey the following equations 



-d v + m 



?¥)m = rn n e^f«(y) 



d y + m^P) f*(y) = m n e^ti(y) (6) 



k 

and the orthogonality relations (taking account of the Z 2 symmetry): 



L dye^f L * m (y)ti(y) = ? dye^ f R * m (y)f*(y) = 5 mn (7) 

L J~L 



where we assume that the length of the orbifold is 2L. 

We solve the above system of coupled differential equations by substituting f^(y) from 
the second in the first equation. Thus we end up with a second order differential equation, 
which can always be brought to a Schrondiger form by a convenient coordinate transfor- 
mation from y to z coordinates related through ^ = e a<yV \ the coordinate transformation 

chosen to eliminate the terms involving first derivatives. Thus we end up with the differ- 
ential equation of the form: 

•j^ + VW*)} /„*(*) = "f /,?(*) (8) 

with t r, M = ^+ i )y)' -^L^ te) (9) 

2[g{z)} 2 2[g(z)] 2 

Here f^{ z ) = fn(v) anc ^ we nave defined v = ^ and g(z) = e a( * y \ The left handed 
wavefunctions are given by []: 

tiiv) = < -^[ d v + m ^f) fniv) (io) 

The form of eq.(8) and (9) are exactly the same as that satisfied by the graviton when 
' r p|] . For any u, we note that before orbifolding the system supports two zero modes 



V ~ 2 



9 Note that it can be shown that the left-handed component obeys also a similar Schrodinger equation 

2[ S (z)P + 2l^)P 



with Vl(z) — V 2\a(zW — ^ 2{aU) ] ' i a " (v) which is the same as Vr f° r v ^ —v 
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(left-handed and right-handed) fl7| . However, the orbifold compactification leaves only a 
chiral right-handed zero mode. 

We note that since the bulk fermion mass, m, is a parameter that appears in the original 
five dimensional Lagrangian its "natural" value is of the order of the five dimensional Planck 
scale M 5 . Now since we assume that k < M 5 (in order to trust our perturbative analysis 
when we consider the configuration gravitationally) it is clear that the "physical" value of 
v is v > 1. However, we will always comment on the behaviour of our results out of this 
region (even for negative values). 

We are particularly interested in the coupling of the bulk spinor to the SM neutrinos 
since this is the way that the neutrino masses will be generated. In order to avoid weak 
scale neutrino masses and lepton number violating interactions we assign lepton number 
L = 1 to the bulk fermion state and thus the only gauge invariant coupling is of the form 

Sy = ~ J d 4 Xy/-g Br {Y 5 L (x)H (x)^ R (x, L Br ) + h.c.} (11) 

where H is the SM Higgs field, L is the SM lepton doublet, H = ioiH* : , g™J is the induced 
metric on the brane and g v i s = det(g™J). The Yukawa parameter Y 5 has mass dimension 
— ~ and thus since it appears a parameter in the five dimensional action, its "natural" value 

i S y, no 1 no 4= 
b 5 V^' 



To obtain canonical normalization for the kinetic terms of the SM neutrino we perform 
the following field rescalings Hq — > e a( - LBr ^H, Lq — > e 3a( - LBr ^ 2 L , where Ls r is the position 
of the brane that SM is confined. This gives 

S Y = - E / d 4 x{y n L(x)H(x)^ n R (x) + h.c.} (12) 

n>0 

where 

y n = e<^-)/ 2 Y 5 fn(L) = {g{z Br )) x " Y 5 fJ(z L ) (13) 

From the above interaction terms we can read off the mass matrix M that appears in 
the Lagrangian as ip u L Niip v R + h.c. where we have defined %l) u L = {yi, . . . , ijj^) and ip v R = 
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("00^) "0?) ■ ■ ■ ) "0n )• The mass matrix for the above class of models has the following form 



M 



( vy vyi 
nil 



V o 








vy n \ 




m n j 



3 Neutrinos in RS model 

For completeness of our analysis, we first briefly review the case of bulk fermion spinors in 
the RS model |17| |. This model consists of one positive (hidden) and one negative tension 



brane (where the SM fields are confined) placed on the fixed points (y = 0, L\) of a S 1 /Z 2 
orbifold (for details see 0). In this case the background geometry is described by eq.(l) 
where a(y) = k\y\. The convenient choice of variable, for the reasons described in the 
previous section, is: 

z^-j— ve[o,£i] (14) 

Since in this model we have (a'(y)) 2 = k 2 and a"(y) = 2kg{z) [S(z) — 5(z — zx)], the po- 
tential appearing in eq.(8) of the Schrodinger equation that the wavefunction of the right- 
handed bulk fermion is (for z > 0): 

V ^ = v -W?-W) 2mz) - S(z - Zl)] (15) 

Here we have defined g(z) = kz + 1 and Z\ = z(Li). 

This potential always gives rise to a (massless) zero mode. It is given by 

ffi') - ^ (16) 

From the above expression it seems that the zero mode is always localized on the positive 
tension brane for all values of v. Nevertheless, by taking the second brane to infinity, we 
find that the zero mode is normalizable in the case that v > \ and that it fails to be 
normalizable when < v < \. The above, as we mentioned, shows that only when v > ~ 
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the zero mode is localized on the first brane. For v — 5 there is no localization and for 
< v < \ it is localized on the negative brane. Another way to see the above is to find 
the coupling of the KK states to mater of a "test" brane as a function of the distance from 
the first (hidden) brane. From eq.(13) we can find that in the case of v > | the coupling 
decreases as we go away from the first brane, on the other hand it is constant when v — \ 
(no localization), and increases when < v < ~ (localization on the second brane). In 
any case as we previously mentioned the "natural" value for v can be considered to be 
greater than unity (having already restricted ourselves in the region v > 0) and thus we 
will assume in the following discussions that the right-handed zero mode is always localized 
on the hidden positive tension brane and we will briefly discuss the rest possibilities. Note 
that the all the following expressions for the masses and the coupling are valid under the 
assumption that v > |, as the results for the rest of the parameter space are different. In 



this case the normalization constant is A ~ y k(y — |). 

Apart from the zero mode we have to consider the left and right-handed KK modes which 
correspond to solutions for m n > 0. The solutions for the right-handed wavefunctions in 
this case are given in terms of Bessel functions 0: 



AJ„ +1 ( ^g(z)) + BJ_ v _i (^g(z) 



(17) 



/ k 

These solutions must obey the following boundary conditions: 

/* W + ^ j?(0) = 

f* '{*!-) + -^r/^i) = (18) 
9{zi) 

The wave functions of the left-handed KK states can be easily extracted from eq.(10). 
The boundary conditions give a 2 x 2 system for A,B which, in order to have a nontrivial so- 
lution, should have vanishing determinant. This gives the quantization of the spectrum. For 

v > I the quantization condition can be approximated by a simpler one: J v _±_ ( mn9 ^ Zl ^ — 

This implies that the KK spectrum of the bulk state is: 

10 Note that in the case that v = N + i, where is an integer, the two linearly independent solutions 
are: Jjv+i and Yjv+i. Although for our calculations we have assumed that v ^ N + h, all the results for 
the mass spectrum and the couplings are valid also in the special cases when v = N + i . 
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m n = Cn k e~ kLl (19) 
(for n > 1), where £ n in the n-th root of J v _i(x) . This means that if one is interested 

in solving the hierarchy in the context of this model, i.e. w = e~ kLl ~ 10~ 15 the mass of 
the first bulk spinor KK state will be of the order of 1 TeV and the spacing between the 
tower will be of the same order. To, summarize the spectrum in this case consists of a 
chiral right-handed massless zero mode and a tower of Dirac KK states with masses that 
start from 1 TeV (if a solution of the hierarchy is required) with ~1 TeV spacing. The 
other important point for the phenomenology is the coupling of the bulk spinors to the SM 
neutrino. It is easy, using eq.(13), to find that the zero mode couples as 



vy Q = vY^k{u- 1 -) (-^-yj a ~ V ^/-l (20) 

since, the hierarchy factor is defined as w = and, as mentioned in the previous section, 
Y 5 ~ ^, v ~ 10 2 GeV and u > \. 

In a similar fashion we can find the couplings of the SM neutrino to bulk KK states. In 
this particular model it turns out that this coupling does not depend on the fermion mass 
or the size of the orbifold and thus it is a constant. By a simple calculation we find that 



vy n 



V2v Y 5 Vk ~\[2v (21) 



Thus from the above we see that the KK tower couples to SM neutrino with a TeV 
strength. In order to find the mass eigenstates and the mixing between the SM neu- 
trino and the sterile bulk modes one has to diagonalize the matrix MM* (actually one 
finds the squares of the mass eigenvalues). By performing the above diagonalization, 
choosing e~ kLl ~ 10 -15 it turns out that the mass of the neutrino will be of the order 
m u ~ 10 2 (lO -15 )*^ (e.g. for v — |, m u ~ 10 -4 eV), and the masses of the bulk states 
are of the order of ITeV with a ITeV spacing. From the last calculations it appears that 
one can easily create a small neutrino mass and at the same time arrange for the desired 
mass hierarchy when v > |. Apart from creating small masses, one has to check that the 
mixing between the SM neutrino and the KK tower is small enough so that there is no 



conflict with phenomenology. It was shown in Ref. |17[ that this can be done for this model 



without fine-tuning. Note that the parameter space: v < ^ (including negative values) is 
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not of interest in the present discussion [^] since it would be impossible to solve the hierarchy 
problem and in the same time to assign small masses to neutrinos. 



4 Neutrinos in " + H — " model 

Since we are interested in studying the characteristics of bulk fermion modes in multi-brane 
configurations we add to the " H — " RS model another positive tension brane where now 
SM fields will be confined. Thus we end up with two different configurations: the " + H — " 

model which will be the subject of this section and the " H h" model which will be the 

subject of the next section. 

The " + H — " model consists of two positive and one negative tension brane. The first 
positive brane is placed on the origin of the orbifold at y = the second (where the SM 
fields are confined), which is freely moving, is place at y — L\ and the negative brane is 
placed at the second fixed point of the orbifold at y = L 2 . 

In the present model the convenient choice of variables is defined as: 



2g«iM— 1 
ki 

e k2(y-Li) + k 1 Ll £ k 1 L 1 _ 
k~ 2 1 IT 



ye [o,Li] 



:22) 



Note the presence of two bulk curvatures, namely h\ and k 2 in this model, which is the 
price that we have to pay in order to place two positive branes next to each other (k± < k 2 
but with k\ ~ k 2 so that we don't introduce another hierarchy. For details see Ref.||33||). In 
terms of the new variables we can find that the potential Vr(z) of the Schrondiger equation 
that corresponds to the present model has the form (for z > 0): 



V R (z) 



f{kl{6{z) - 6{z - z x )) + k 2 2 (6(z - Zl ) - 6(z - z 2 ))) 



hSiz) 



(fc 2 -fci) 

2 



5(z - zi) - k 2 5(z - z 2 ) 



(23) 



since cr"(y) = 2g(z) k\5(z) + 



(fc 2 -fci). 



S(z — Zi) — k 2 5(z — z 2 ) and (a'(y)) 2 = k\ for y G 



11 This could be of particular interest if one uses the above mechanism to localize SM fermions on the 
negative tension brane and in the same time solving the hierarchy problem (e.g. see Rcf.[^l|). 
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[0, L\] and (a'(y)) 2 = k\ for y G [Li, L 2 \. The function g(z) is defined as: 



9{z) 



k x z + 1 



z e [0,^] 



/c 2 (-2 - ^l) + fci^i + l z e [z b Z 2 ] 



(24) 



where z$ = 0, zi = z(£i) and z 2 = z(L 2 ) are the positions of the branes in terms of the 
new variables. 

This potential always gives rise to a (massless) zero mode whose wavefunction is given 

by 

A 



/<?(*) 



(25) 



The discussion of the previous section about the state localization applies in this model 
as well . For v > \ the zero mode is localized on the first brane. For the case v — \ there is 
no localization again. For v < \ it is localized on the negative tension brane, as expected. 



In the case v > \ we find that the normalization factor of the zero mode is A ~ a / k\ 



which is the same as in the case of RS for k = k\ (not surprisingly since it is strongly 
localized on the first brane). 

The wavefunctions for the right-handed KK modes are given in terms of Bessel functions. 
For y lying in the regions A = [0, L{\ and B = [L 1: L 2 ], we have: 



x£/(«) 



A 
B 



£(£) 
fci 

k 2 



AoJ 



+ B 2 J_ V _, 



(26) 



with boundary conditions: 



R 1 1 



i? w + ^fm = o 

f«(zi + ) - fj(zi-) = 
f*(zi) = 



2v (k 2 -k x \ ~ eR 



9(zi) 



hp 



/„* \z 2 -) + ^mzi) = o 



(27) 
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The above boundary conditions result to a 4 x 4 homogeneous system for A 1 , Bi, A 2 
and B 2 which, in order to have a nontrivial solution, should have a vanishing determinant. 
This imposes a quantization condition from which we are able to extract the mass spectrum 
of the bulk spinor. The spectrum consists, apart from the chiral (right-handed) zero mode 
(massless) which was mentioned earlier, of a tower of Dirac KK modes. 

In this case in order to provide a solution to the hierarchy problem we have to arrange 
the distance between the first two branes so that we create the desired hierarchy w. In the 
present model we have an additional parameter which is the distance between the second 
and the third brane x = k 2 {L 2 — Li) For the region where x >1 we can find analytically 
that all the masses the KK tower (including the first's) scale the same way as we vary the 
length of the orbifold L 2 : 



where ( n is the n-th root of J u _i(x) = 0. 

In the region x < 1 the previous relation for the mass spectrum breaks down. This 
is expected since for x = (L 2 = Li) the " + H — " model becomes " H — " (RS) and the 



RS condition. So for < x < 1 the quantization condition (and thus the mass spectrum) 
interpolates between the previous two relations. 

Let us now turn to the coupling between the SM neutrino which lives on the second 
positive brane with the bulk right-handed zero mode and the rest of the KK tower. We can 
easily derive that zero mode couples in that same way as in the RS case (the normalization 
of the zero mode is approximately the same): 



m n = Cn wk 2 e 



-k 2 L 2 



(28) 





(29) 



where, as we previously mentioned, Y 5 ~ -|=, v ~ 10 2 GeV, v > \ and w = -At- On the 



other hand the coupling of the SM neutrino to bulk KK states is given by: 




(30) 



where ( n is the n-th root of J u _i(x) = 0. 
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All approximations become better away from v — ~, x — 0, and for higher KK levels. 
Note the strong suppression in the coupling scaling law. This rapid decrease, which is 
distinct among the models that we will consider, also appears in the coupling (to matter) 
behaviour of the graviton KK states and a detailed explanation can be found in Ref. ||33||. 

Thus from the above we conclude that for v > | the phenomenology of this model 
resembles, in the general characteristics, the one of RS. Of course in the present model 
there is an extra parameter, x, which controls the details of the masses and couplings of the 
KK states. Since the zero mode coupling is independent of x the general arguments of the 
previous section about creating small neutrino masses apply here as well, at least for small 
x . By increasing x we make the KK tower lighter, as we see from eq.(28), but we avoid 
large mixings between the SM neutrino and the left-handed bulk states due to the fact that 
the coupling between the SM neutrino and the right-handed bulk states drops much faster 
according with eq.(30). 

Note that in the case v < ~ (negative values included) the zero mode will be localized 
on the negative brane and thus one could arrange the parameter x so that the exponential 
suppression of the bulk fermion's zero mode coupling on the second brane is such that gives 
small neutrino masses. Thus in this case it seems that we are able to solve the hierarchy 
problem by localizing the graviton on the first positive brane and in the same time create 
small neutrino masses by localizing the bulk fermion zero mode on the negative brane. 
Nevertheless, one should make sure that no large mixings are induced in this case. 



5 Neutrinos in " H \-" model 



We now turn to examine bulk spinors in the " H \-" model, which was analyzed in detail 



in Ref. p9|, pz| . The model consists of two positive tension branes placed at the orbifold 
fixed points and a third, negative brane which is freely moving in-between. SM field are 
considered to be confined on the second positive brane. Of course the presence of a moving 
negative brane is problematic since it gives rise to a radion field with negative kinetic term 
(ghost state) Q in the gravitational sector. Nevertheless we are interested in the 
general characteristics of this model . The interesting feature of this model is the bounce 
form of the warp factor which gives rise to an ultralight graviton KK state as described in 
Ref. . It was shown in Ref. [ 3~3] that exactly this feature, of a bounce in the warp factor, 
can be reproduced even in the absence of negative branes in the " + +" model where this 
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is done by sacrificing the flatness of the branes (the spacetime on the branes in this case 

is AdS 4 ). Thus we will handle the " H V" as a toy model since in this case there can be 

simple analytical calculations of the coupling etc. The general characteristics will persist 
in the " + +" case. 

We are interested to see if this configuration as well as an ultralight graviton supports 
an ultralight spinor field. In order to see this we should check the form of the potential of 
the Schrondiger equation that the right-handed component obeys. We can easily find that 
the potential is (for z > 0): 

Vr{z) = V{ 2\g{z 1 )? ~ ^¥) 2k m + 6(Z " Z2) " " Zl)] (31) 

since (a'(y)) 2 = k 2 and <r"(y) = 2kg(z) [S(z) + S(z — z 2 ) — S(z — Zi)]. The convenient choice 
of variables in this case is: 



z = < 



2e kL l —e 2kL l~ ky — 1 
k 



y e [Li,L 2 ] 
^ ye[0,Li] 



(32) 



and the function g(z) is defined as g(z) = k \z\ — \\z\ — zi|} + 1, where z\ = z{L\). 

As in the previous cases, the above potential always supports a (massless) zero mode 
with wavefunction of the form: 

f ° {z) ^ WW (33) 

In this case the different localization behaviour as a function of v is the following: For 
v > | the zero mode is localized on the positive branes (thus fails to be normalizable when 
we send the right positive brane to infinity but is normalizable when we send both negative 
and positive to infinity). For the case v < \ the localization of the zero mode is on the 
negative tension brane, as expected. In the case v > | and for strong hierarchy w we find 

that the normalization factor of the zero mode is A ~ \Jki{v — |) (Note that in the case of 
"weak" hierarchy one should be careful with the assumptions on which the approximations 
are based on e.g. for w — 1 the result must be divided y/2). For the KK modes the solution 
is given in terms of Bessel functions. For y lying in the regions A = [0, Li] and B = L 2 ], 
we have: 
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(34) 



with boundary conditions: 



fj '(0 + ) + JffiffW = 

/:v)-/rM-|^/^i) = o 

w) - -^%&) = ° ( 35 ) 

The boundary conditions give a 4 x 4 linear homogeneous system for A\, B\, A 2 and 
B 2 , which, in order to have a nontrivial solution should have vanishing determinant. This 
imposes a quantization condition from which we are able to extract the mass spectrum of 
the bulk spinor. The spectrum consists, apart from the chiral (right-handed) zero mode 
(massless) which was mentioned earlier, by a tower of Dirac KK modes. Nevertheless due 
to the fact that there are two positive tension branes present in the model there are now 
two "bound" states in a similar fashion with Ref.pUL (for v > |). One is the the right- 



handed zero mode which is massless and it is localized on the positive brane placed at the 
origin of the orbifold and the second is the ultralight right-handed first KK state which is 
localized on the second positive brane placed at the other orbifold fixed point. This can be 
seen by examining the mass spectrum and the coupling behaviour of the first KK state in 
comparison with the rest of the tower. 

Firstly let us examine the mass spectrum. In the case that we have a hierarchy w (where 
w = -^-^y = e~ a ( L ^) we can find appropriate analytical expressions for the mass spectrum. 

For the first KK state 

mi = VAu 2 - 1 kw e - {v+ ^ )x (36) 

and for the rest of the tower 

m n+ i = £ n kw e~ x n = 1, 2, 3, . . . (37) 
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where £22+1 is the (i+l)-th root of J u _i(x) (i = 0, 1, 2, . . .) and £ 2 i is the i-th root of J u+ i (x) 

(i = 1,2,3, .. .). The above approximations become better away from the v = \ , x = 
and for higher KK levels n. The first mass is manifestly singled out from the rest of the 
KK tower as it has an extra exponential suppression that depends on the mass of the bulk 
fermion. By contrast the rest of the KK tower has only a very small dependence on the mass 
of the bulk fermion thought the root of the Bessel function £ n = £ n {v) which turns out to 
be just a linear dependence in v. Note there is a difference between the graviton ultralight 
state (discussed in |29|, |3"2|| ) and this spinor state: In the case of gravity the unltralight KK 
state the mass scales as a function of x was e~ 2x , on the other hand the scaling law in the 

case of the ultralight spinor is of the form e~( u+ ^ x . From the above it seems that the latter 
can be done much lighter that the graviton first KK state for a given x by increasing the 
parameter v. This is easy to understand since the role of the mass term, with the kink or 
multi-kink profile, is to localize the wavefunction f(z). By increasing the parameter v all we 

do is to force the absolute value of the wavefunction of the first KK state and the massless 
right-handed zero mode to become increasingly similar to each other: For example, in the 
symmetric configuration, the the difference between the zero mode and the first KK state 

wavefunctions comes from the central region of the " H V" configuration, where the first 

KK state wavefunction is zero (since it is antisymmetric) thought the zero mode's is very 
small due to the exponential suppression of the wavefunction, but non zero. By increasing 
v we force the value of the zero mode wavefunction at the middle point to get closer to zero 
and thus to resemble even more the first KK state, something that appears in the mass 
spectrum as the fact that the mass of the first KK state is approaching to zero. On the 
other hand the mass eigenvalues that correspond to the rest of the tower of KK states will 
increase linearly their mass by increasing the v parameter since those are not bound states 
(the first mode has also such a linear dependence in v but it is negligible compared with 
the exponential suppression associated with v ). 

Now let us turn to the behaviour of the coupling of the zero mode and the KK states to 
matter living on the third (positive) brane. As in the previous cases the right-handed zero 
mode couples to SM left-handed neutrino as 

v y = v Y 5 \Jk(v-~) (-^j 2 ~ v y / ^| w v ~* (38) 

since Y$ ~ -7=. From the above relationship we see that the coupling of the zero mode 
to SM neutrino will generally be suppressed by the hierarchy factor to some power, the 
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power depending on the bulk fermion mass. This way one may readily obtain a very small 
coupling. The coupling of the zero mode is independent of x. This is another way to see 
the localization of this mode on the first brane (the normalization of the wavefunction is 
effectively independent of x). Since this model supports a second "bound state" (first KK 
state) which is localized on the second brane, we expect something similar to occur in the 



coupling behaviour of this state. Indeed, similarly to the graviton case |29], |32|, we can show 
that the coupling of this state to the SM neutrino for fixed w is constant, i.e. independent 
of the x parameter. Taking in account the result of the graviton KK state a\ = w l Ipi and 

by comparing the graviton-matter and spinor matter coupling we can easily see that the 
coupling of this special mode will be of the order of the electroweak scale: 



vyi - y u ~\ v ( 39 ) 

Let us now consider the coupling of the rest right-handed KK states to the SM neutrino. 
We find that 



v y n ~ Jv- -v e x (40) 

for n = 0, 1, 2.... From the above relationship we see that the rest of KK states will generally 
have exponentially suppressed coupling compared to the first special state. 

The appearance of this special first ultralight and generally strongly coupled KK state, 
as in the graviton case, is going to have radical implication to the phenomenology of the 
model. Let us consider the following example suppose that v = | and that we also require a 
hierarchy of the order: w ~ 10 -15 . In this case the zero mode's coupling is vyo ~ id0~ 15 ~ 
1CT 4 eV a result independent of the x parameter. On the other hand one can check that 
the rest of KK tower will have masses m n ~ 10 3 e~ x GeV (for n = 2,3...) with coupling 
vy n ~ ve~ x . Up to this point the phenomenology associated with this model is similar to the 
RS case i.e. tiny coupling of the right-handed and generally heavy KK states with relatively 
strong coupling. However, taking in account the special KK state, we have the possibility 
of obtaining a much lighter state with large coupling (effectively independent of how light 
this state is). Having a light sterile state whose right-handed mode has strong coupling 
to the SM neutrino is potentially dangerous. In such a case we find that the dominant 
contribution to the mass eigenstate of the lightest mode (neutrino) v p h ys will come from 
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the left-handed component of this special sterile mode and not the weak eigenstate v. Of 
course something like this is not acceptable since there are strict constrains for the mixing 
of SM neutrino to sterile states. Since the mass spectrum depends exponentially on the 
the x parameter which determines the distance between the branes, the above argument 
impose strong constraints on it's possible values. 

Finally, in the case that v < | the bulk right-handed zero mode is localized on the 
negative tension brane. In this case a new possibility arises: By localizing the graviton 
wavefunction on the first brane we can explain the SM gauge hierarchy (by setting w to 
the desired value ) and by localizing the bulk fermion zero mode on the negative tension 
brane to induce small neutrino mass (for appropriate value of the x parameter) for the SM 
neutrino which is confined on the right positive tension brane. Note that for v < | there 
is no special bulk spinor KK state and thus there is no immediate danger of inducing large 
neutrino mixing from such a state. However, the presence of the ultralight graviton KK 
state is restricting our parameter space as following: In order to solve the gauge hierarchy 
problem (assuming the SM on the third brane) we have to fix the one parameter of the 
model: w ~ 10~ 15 . Since the fermion zero mode is localized on the intermediate negative 
brane we have to arrange the distance between this and the third brane (for given u), x, so 
that the coupling is sufficiently suppressed in order to give reasonable neutrino masses. This 

implies that e~\ u ~^ x ~ 10~ 13 , if one considers the mass of the neutrino of the order of 10 _1 
eV. From the bounds derived in |!2] we find, for k ~ 10 17 GeV, that in order the ultralight 



graviton KK state not to induce modifications of gravity at distances where Cavendish 
experiments take place and not to give visible resonances to e + e~ — > fi + processes, we 
should have 4.5 < x < 15 or x < 1. The latter implies certain restrictions to the values of v: 
1.5 < —v < 6.2 or —v > 29.4. In the above regions it is possible to simultaneously create 
the gauge hierarchy and small neutrino masses consistently, with the mechanism described 
earlier. 

6 Neutrinos in " + +" model 

As we mentioned in the previous section the " + +" model mimics the interesting char- 
acteristics of the " H V" model without having any negative tension brane. Thus since 

the warp factor has a bounce form this model also supports an ultralight graviton as was 



shown in Ref.|34|, [35], [36], [57] . According to the previous discussion, we should expect that 



the model will support a ultralight sterile neutrino as well. This can be easily shown again 
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by considering the form of the potential of the differential equation that the right-handed 
component is obeying, which will again turn out to be of the same form as the graviton. 

Now it turns out that for the construction of such a " + +" configuration, it is essential 
to have AdS± geometry on both branes (for details see Ref.|34j|). Thus in this case the 
background geometry is described by: 

e -Mv) 

ds 2 = H2x 2,r > Vnvdx> J 'dx v - dy 2 (41) 



1 - 



4 / 

where the corresponding inverse vielbein is given by 

H 2 2 rr2 2 tt2 2 rx2 2 

Eg = diag(e°W(l - - ^-),e*M(l - ^-),e"<«>(l - ^-), 1). (42) 

Since now the brane is no longer flat the previous calculations for the action will be 
slightly modified. We briefly discuss these modifications. Following the same steps of the 
flat case, we write = \I/ R + where ^ r> l = |(1 ± 7s)^- Since the connection part of 
the Lagrangian again doesn't give any contribution (since the vielbein is again symmetric) 
the action becomes: 



S = j^dd^x dy{e~ 3(7 E£ {^ Lil a d A ^ L + ^ R ij a d A ^ R ) - e^m^- (f L * fl + # fl tt L ) 

-~ [y L (e- 4 °d y + d y e- 4 ^R ~ ^R{e- Aa d y + d y e~^)m L ] }(43) 

where E£ = (1 — ^-f-)5a with (a,A=0, 1,2,3) is the induced vielbein and G the determinant 
of the induced metric. For convenience and in order to be able to use results of RefJB4 



we set A(y) = e °w where 

= cosh(t(y - |y|)) 
cosh(% ) 

is the equivalent "warp" factor in this case, which is found by considering the configuration 
gravitationally. From the above relation it is clear that the "warp" factor has the desired 
bounce form, with a minimum at y$. The position of the minimum, something that it is 
going to be important for the phenomenology of the model, is defined from the relationship 
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tanh(/cy ) = pqS where V\ is the tension of the first brane and A is the five dimensional 
cosmological constant. As we mentioned in the introduction, the profile of the mass term, 
^^y, in the present model is not a simple combination of ^-functions. In particular there 

are two ^-function profiles near the orbifold fixed points which give rise to the positive 
branes, but there is also an intermediate kink profile of the form — tanh(/c(y — yo)) which 
is associated with the presence of the bounce (this could give rise to a "— " brane as a limit 

, resulting to the familiar " H \-" configuration). Note that even though there is no brane 

at the position of the minimum of the "warp" factor the kink profile is expected to act in 
the same way, and thus induce localization of the fermion zero mode in specific regions of 
the parameter space, exactly as in the " H h" model. 

As in the flat-brane case, we can decompose the left-handed and right-handed fermion 
fields into KK states with nontrivial profile wavefunctions f%,fn (in respect to the fifth 
dimension) in order to be able to bring the Lagrangian into the form 

S = J2 d 4 x^§{E^tp n (x)i-f a d A ipn(x) - m n ^ n (x)il) n (x)} (45) 

n J 

where the wavefunctions fn(y), fniv) should obey the following equations 



(-^ + j^)ti(y) = rn n A-\y)f«(y) 
(d y + f«(y) = m n A-\y)f^y) (46) 

with the following orthogonality relations: 

\ L _ L dyA-\y)j L m \y)] L n {y) = J^dyA-^f^f^y) = 5 mn (47) 

Again we solve the above system of differential equations by finding the second order 
differential equation that it implies for the right-handed component of the spinor. It is 
always possible to make the coordinate transformation from y coordinates to z coordinates 
related through: ^ = A~ l {y) and bring the differential equations in the familiar form: 
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(48) 



where we have defined f^{z) = f^(y) 



with Vr(z) 



^A(y)A"(y) + ^(A\y)f 
u 2 k 2 u(u+l)k 2 



COS' 



(i 



z - z 



— kv 



tanh(% )5(,) + ^(^-yo)) cosh(fc(L-y )) _ " 

cosh (%o) 



(49) 



with k defined as k = — ry; — r. The new variable z is related to the old one y through the 

cosn(feyo) 

relationship: 



z = sgn(y) 7 



arctan (tanh( — — ) J + arctan (tanh(-^ 



(50) 



Thus in terms of the new coordinates, the branes are placed at z± — and z L , with the min- 
imum of the potential at z — | arctan (tanh(^p)) . Also note that with this transformation 

the point y = oo is mapped to the finite point Zoo — | | + arctan (tanh(^p 



We can now proceed to the solution of the above equations for the right-handed compo- 
nents, while the left-handed wavefunctions can be easily evaluated using eq.(10) (taking in 
account the definition A = e~ CT( ^ and the change of variables). The zero mode wavefunction 
is given by: 

C 



/«?(*) 



(51) 



[cos(k(z - \z\))] v 

where C is the normalization factor. If we send one of the two branes to infinity (i.e. 
z\ — > Zoo) an d at the same time keep zq fixed we find that the zero mode is normalizable 
only in the cases where v < |. In the other cases [y > |) the wavefunction fails to be 
normalizable due to the fact that it is too singular at z^. Note though that in the case 
where v < | the first KK will not be special (i.e. will have almost the same behaviour as 
the rest of the KK tower). We also note that again the zero mode is chiral i.e. there is 
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no solution for f£ when m = that can satisfy the boundary conditions (antisymmetric 
wavefunction). From the above we see that the localization behaviour of zero mode in the 

" + +" model is the same as in the " H y" model. Note that for v < | the zero mode will 

be localized near y despite the absence of any brane at that point. This is because, as we 
mentioned, the — f ac tor, which can be considered as the vacuum expectation value of 

a scalar field, has a kink profile in the neighbourhood of yo which induces the localization 
(this kink becomes the negative tension brane in the flat brane limit). By considering cases 
with m n 7^ 0, we find the wavefunctions for the KK tower : 



f*(z) = cos" +1 (fc(|z| - zo)) [Ci F(a n: b n , ±; sm 2 {k{\z\ - z ))) 

+ C 2 \ sm(k(\z\ -zq))\ F(a n + ±,b n + I, I; sin 2 (k(\z\ - z ))) 



(52) 



where 



2 2 V V k 



+ * (53) 



2 V V A; 

The boundary conditions are given by: 

(0 + )+^tanh(^o)/,f(0) = 

the above conditions determine the mass spectrum of the KK states. By studying the mass 
spectrum of the KK states it turns out that it has a special first mode similar to the one 

of the " H h" model as expected. For example, for the symmetric configuration (w — 1), 

by approximation we can analytically find the following expressions for the mass of this 
special state : 

mi = 2V4z/ 2 + 3 k (e- kyo Y + * (55) 
In contrast, the masses of the next levels are given by the formulae: For odd states 



m n = 2yJ(n + l)(n + l + u) k e~ kyo (56) 
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with n = 0, 1, 2, and for even states 



m n = 2^ (n + |) (n + | + i/) A; e"*" (57) 

with n = 0, 1,2, .... 

Again the mass of the first KK state is manifestly singled out from the rest of the KK 
tower as it has an extra exponential suppression that depends on the mass of the bulk 
fermion. The above characteristics persist in the more physically interesting asymmetric 

case [w « 1). In this case we find mi ~ m\{v) kw (e -fcr ) l/+ 5 and m n ~ m°(i/) kw e~ kx 
where x = L\ — yo is the distance of the second brane from the minimum of the warp factor, 
and rrii ~ 1j m n * s linear to n (for large values) and has a ~ y/v dependence. The first 

KK state is localized on the second positive brane and, as in the case of " H h" model, 

its coupling to SM neutrinos remains constant if we keep the hierarchy parameter w fixed. 

The phenomenology of this model will be similar to the " H V" model and thus we do not 

consider it separately. 



7 Bigravity and Bulk spinors 



Bigravity [^, |32], |34| (multigravity |30], |33fl ) is the possibility that gravitational interactions 



do not exclusively come from a massless graviton , but instead they can be the net effect of 
a massless graviton and one or more KK states (continuum of KK states) or even a single 
massive KK mode [[34], |35| without conflict with General Relativity predictions fl40 , 4l| . This 



is based on the different scaling laws of the mass between the first and the rest of KK states. 
Since the first graviton KK state has mass with an additional exponential suppression we 
can realize the scenario that the first KK state is so light that its wavelength is of the order 
of the observable universe and thus any observable effect of its non-vanishing mass to be out 
of the experimental reach, and in the same time the rest of the KK tower has masses above 
the scale that Cavendish experiments have tested Newtonian gravity at small distances. 

In the two previous sections we have shown that in the case of models, where Bigravity 
can be realized, there is also an ultralight KK that corresponds to the bulk spinor assuming 
the existence of the mass term for the bulk spinor that appears in eq. (2) (with v > i) . The 
subject of this section is to investigate if the two above possibilities are compatible: Can 
we have a Bigravity scenario and a consistent neutrino phenomenology? For the sake of 
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Figure 2: On the left, the mass spectrum of the graviton (G) (first column) and bulk 

spinor (S) (second column) KK states in the " H h"-Bigravity model. On the right, for 

comparison, the corresponding spectrum for the case of RS model. The figures are not in 
scale and the details of the bulk spinor mass spectrum depend on the additional parameter 
v and thus in the above figure we have assumed v ~ | (for higher value of v we could have 
the fermion first KK state to be lighter that the first graviton KK state). 



simplicity the discussion below will be concentrated to the " H V" model but, analogous 

arguments should apply to the case of " + +" model. 

Let us review briefly the Bigravity scenario. The graviton ultralight first KK state has 



mass = 2wk e 2x and coupling a^' 



A . The rest of the KK tower has masses : 

wMpi 



rn^+i = £n kw e~ 



n 



1,2,3, 



(58) 



where (,2i+i is the (i + l)-th root of Ji(x) (i = 0, 1, 2, . . .) and ^ is the i-th root of Jz{x) 
(i = 1,2,3, . . .). The couplings of these states scale as a n oc e~ x . In order to achieve the 
Bigravity scenario, as we mentioned in the beginning of this section, we have the following 

constrains on the range of masses of the KK states: mf^ < l(T 31 eV or mf ] > l(T 4 eV 
(where Planck suppression is considered for the "continuum " of states above 10~ 4 eV). Our 

exotic scheme corresponds to the choice ~ 10~ 31 eV and > 10~ 4 eV. In this case, 
for length scales less than 10 26 cm gravity is generated by the exchange of both the massless 
graviton and the first KK mode. This implies, (taking into account the different coupling 
suppressions of the massless graviton and the first KK state) that the gravitational coupling 
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as we measure it is related to the parameters of our model by: 



1 




) 



1 



2 => Mpi « wM 5 



(59) 



(»M 5 ) 



We see that the mass scale on our brane, wM 5l is now the Planck scale so, although the 
"warp" factor, w, may still be small (i.e. the fundamental scale M 5 >> Mpi anc k), we do 
not now solve the Planck hierarchy problem. Using the equations for the mass spectrum 
and assuming as before that k ~ M 5 , we find that mi = 2 kw e~ 2x ~ Mpi anc ^e~ 2x . For 
777,1 = lCT 31 eV we have m 2 « 10~ 2 eV. This comfortably satisfies the bound m > 10 _4 eV. 

Now let us see what this implies for the neutrino physics. Let us first consider the case 
where v > |. In this mentioned above, there will be also an ultralight bulk fermion 

KK mode. By forcing the graviton first mode to have a tiny mass we also force the first 
spinor KK state to become very light (even in the best case where v — > ~ this state will 
have mass of the order of a fraction of eV. For larger v becomes even lighter) since they are 

related through: m\ ~ e^~^ x mf. This is unacceptable since this mode has a constant 
coupling of the order of the weak scale, v , which will induce large mixing of the neutrino 
with the left-handed component of this state. However note that in the limit v — > ^ the 
special first fermion KK mode will become a normal one losing it's localization and thus 
the above may not apply. Unfortunately, this is not the case as in this limit the fermion 
zero mode is delocalized (for v = | the coupling is constant across the extra dimension) 
giving no possibility of inducing small neutrino masses (since the compactification volume 
is very small)Q. 

Despite the severe constraints in the above scenario, due to the presence of the ultralight 
bulk fermion KK state, by no means the Bigravity scenario is excluded in the the case of 
v > g since one can always consider the possibility of placing the SM on a brane (with 
tiny tension so that the background is not altered) between the negative and the second 
positive brane, so that the coupling of the bulk fermion first KK state to SM neutrinos is 
sufficiently small while a part of gravitational interactions will still be generated from the 
ultralight graviton first KK state. 

Let us turn now to the case < v < ^. In this case there is no special bulk spinor KK 
state and thus the above arguments do not apply. In this case the bulk fermion zero mode 

12 Trying to use this window of the parameter space seems like a fine-tuning though since one needs to 
dclocalize the fermions first KK state enough in order to have small coupling to SM neutrino and on the 
other hand to prevent the derealization of the fermion zero mode with the same mechanism (y — > i). 
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is localized on the negative tension brane. Nevertheless, in order for the Bigravity scenario 
to be possible we should have x ~ 60. If we now try to generate neutrino masses with the 
mechanism described in the first section we will find that the coupling of the zero mode 
to SM neutrino will be by far too small to provide consistent results. Thus in this case 
although Bigravity is realized we cannot generate neutrino masses consistently. 

The case of v = | is of no interest since in this case there is no localization (the fermion 
zero mode has constant coupling across the extra dimension). The case where v < 
resembles the < v < ~ case, with the only difference that the localization on negative 
branes will be even sharpner, making the situation even worse. 



8 Discussion and conclusions 

We have studied bulk fermion fields in various multi-brane models with localized gravity. 
The chiral zero model that these models support can be identified as a right-handed sterile 
neutrino. In this case small neutrino Dirac masses can naturally appear due to an analogous 
(to graviton) localization of the bulk fermion zero mode wavefunction without invoking a 
see-saw mechanism. For models in which the localization of the fermion zero mode is 
induced by the same scalar field that forms the branes the localization behaviour of this 
mode can resemble the graviton's at least in a region of parameter space. The latter implies 

that the " H h", " + +" models can support, in addition to the ultralight graviton KK 

state, an ultralight localized and strongly coupled bulk fermion KK mode. This fermion 

state, when exists, imposes even more severe constrains on the parameter space of " H h" 

and " + +" models. In the case that one requires the Bigravity be realized the light fermion 
KK mode can induce too large mixing between the neutrino and the KK tower and thus it 
restricts even more the allowed parameter space of the relevant models. 



As a general remark we see that the appearance of Multi-localization [[42] in the Multli- 
Brane world picture and its relation to the existence of ultralight states in the KK spectrum 
is not a characteristic of the graviton only, but can also occur in spin 0, | and | bulk states 
with appropriate bulk mass terms Vector fields on the other hand can have also 

ultralight states but in this case a coupling to a bulk scalar field is necessary in order to 



achieve localization [23 



The appearance of ultralight KK states with the additional characteristic, in the case 
of spin i fields, where the left-handed and the right-handed components are localized in 
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different places in the extra dimension, can give new interesting possibilities which will be 
discussed in another publication j£| . 
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